Naked Singularity Formation In Brans-Dicke Theory 



Amir Hadi Ziaic/'B Khcdmat Atazadeh^-Q and Yascr Tavakoli 2, [l 

'Department of Physics, Shahid Beheshti University, Evin, Tehran 19839, Iran 
2 Departarnento de Fisica, Universidade da Beira Interior, 
Rua Marques d'Avila e Bolama, 6200 Covilha, Portugal 
(Dated: May 26, 2010) 

Gravitational collapse of the Brans-Dicke scalar field with non-zero potential in the presence of 
matter fluid obeying the barotropic equation of state, p — wp is studied. Utilizing the concept of 
expansion parameter, it is seen that the cosmic censorship conjecture may be violated for w = — | 
and w = — | which correspond to cosmic string and domain wall, respectively. We show that 
physically, it is the rate of collapse and the presence of Brans-Dicke scalar field that govern the 
formation of a black hole or a naked singularity as the final fate of dynamical evolution and only 
for these two cases the singularity can be naked as the collapse end state. Also the weak energy 
' condition is satisfied by the collapsing configuration. 
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' I. INTRODUCTION 

<N . 

The process of gravitational collapse has been studied since about 60 years ago as a solution of Einstein field 
equations. When a sufficiently massive star many times the size of the Sun exhausts its nuclear fuel without reaching 
' an equilibrium state such as a neutron star or white dwarf, it collapses under its own pull of gravity at the end of 
5-h , its life cycle. Therefore, gravity overtakes and dominates the other three forces of nature, in particular, the weak 
i and strong nuclear forces, which generically provide the outward pressure in a star to balance it against the inward 

pull of gravity. In such ultra-strong gravity regions the densities and spacctime curvature diverge and a spacetime 
C*~) ■ singularity will be born, either hidden within an event horizon (a black hole) or visible to the external universe (a 
naked singularity), as predicted by the singulary theorems in general relativity. The visibility or otherwise of the 
singularity to outside observers is determined by the causal structure of the dynamically developing collapsing cloud 
as governed by the Einstein's field equations. When the internal dynamics of the collapse delays the formation of the 
horizon, the singularity becomes visible, and may communicate physical effects to the external universe 0]. 

Up until now a great deal of effort has gone into the study of the nature of singularities and large classes of the 
solutions of Einstein's field equations treating singularities have been represented. The general and exact class of 
solutions of Einstein's field equations describing spherically symmetric pressureless matter (dust) for motion with no 
particle layers intersecting, independent of the homogeneity assumption, was originally introduced by Lemaitre^ 
in an attempt to describe cosmology, which was further developed and studied by Tolman 0, [H an( i Bondi 
This class could be used to model the gravitational collapse of matter from general inhomogeneous initial conditions 
and the end state of gravitational collapse of a massive star can be studied within this framework Q. Contrary 
to the collapsing Fricdmann case, in which the physical singularity occurs at a constant epoch of time, namely at 
t = 0, the singular epoch (the time t = to(r) at which the area radius of the collapsing shell of matter at a constant 
value of the co-moving coordinate r vanishes) in LTB model is a function of r as a result of inhomogeneity in the 
matter distribution. This model can exhibit two kinds of naked singularity: a shell-crossing 0, Q and a shell- focussing 
smg ularity 0,1 El. These kinds of naked singularity arc generic for spherical dust collapse and with proper choice 
of initial data they are globally naked. The main reason which causes this model to be deficient is that it neglects 
the pressure, which is likely to diverge at these singularities (where the density diverges) 0. A special case of these 
classes of solutions which has served as the basic paradigm in black hole physics is the Oppenhcimcr-Synder study of 
a completely homogeneous, pressure-free and spherically symmetric dust cloud collapse, where a dust cloud undergoes 
a continued collapse which commences from regular initial data (i.e., there is no trapped surfaces forming at the initial 
spacelike surface from which the collapse begins and the light rays can escape from the surface of the star to faraway 
observers) to form a black hole [1, [ll|, [l2j ■ Most of our knowledge about the gravitational collapse is still based on 
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that model. Also it describes well the formation of the horizon and evolution of the central, space-like singularity. 
Oppenheimer and Snyder analyzed the causal structure of the solution by considering in particular an observer on the 
surface of the dust cloud sending signals to a faraway stationary observer at regularly spaced intervals as measured 
by his own clock. They discovered that as the radius of the dust cloud approaches 2M, the spacing between the 
arrival times of these signals to the faraway observer becomes progressively longer, tending to infinity. This effect 
has since been called the infinite redshift effect. The observer on the surface of the dust cloud may keep sending 
signals after its radius has become less than 2M, but these signals can not escape from this region since the speed of 
light is the limiting propagation velocity for physical signals, within a finite affine parameter interval, these signals 
reach a true singularity at r = [lj| . This singularity is formed inside a black hole, as a final state of the collapse 
process and also accords to the concept of the cosmic censorship. The cosmic censorship conjecture which is initially 
represented by Penrose [hij , explain the properties of the final singularity of a gravitational collapse. The main result 
of this conjecture is that, all spacetime singularities arising from regular initial data that appear in a gravitational 
collapse (in an asymptotically flat universe) are always surrounded by an event horizon and hence invisible to outside 
observers (no naked singularities). Moreover in the strong version of this conjecture, such singularities are not even 
locally naked, i.e. there is not any timelike or null geodesies which can emerge from these singularities and go to future 
infinity Q • This hypothesis plays a fundamental role in both the theory and applications of black hole physics and 
has been recognized as one of the most important open problems in classical general theory of relativity. Up to now 
many exact solutions of Einstein's field equations with several kinds of field-sources which admit naked singularities 
depending upon the nature of the initial data and the kinematic al p roperties near the singularity have been considered. 
The models studied so far include the collapse of scalar fields [15|, [l6[ as well as other matter sources including dust 
SlSIIi, radiation [13, [3 El, perfect fluid 0, [IE [13] , imperfect fluids [13, S3, and null strange quark fluids £3,^. 

One would like to inquire about what are the possible physical factors operating during the process of a continual 
collapse of a massive matter cloud which results in the formation of a naked singularity or a black hole as the 
collapse end state. Such an investigation should help us understand much better the physics of black holes and naked 
singularity formation in gravitational collapse. It was recently shown [2^, [24[ that for the spherical dust collapse the 
shearing effects and inhomogeneity present within the collapsing cloud do play a crucial role in delaying the formation 
of the trapped surfaces and the apparent horizon. These effects could in fact make the geometry of apparent and 
event horizons distorted sufficiently, which exposes the singularity to the external observers [23j . 

General relativity is not the only gravitational theory which explain the gravitational phenomena. There are 
alternative theories of gravity that can explain the gravitational phenomena and lead to many interesting results 
and physical interpretations. One of those, is the scalar-tensor theory of gravity, which has been actively studied as 
an alternative and successful theory of gravity. Recently, they have attracted much attention, in part because they 
emerge naturally as the low-energy limit of many theories of quantum gravity, such as Kaluza-Klein theories [25i [2q | 
and supersymmetric string theories [25|, [27j ■ These theories are also important for "extended" cosmological inflation 
models[25|,[28|, in which the scalar field allows the inflationary epoch to end via bubble nucleation without the need for 
fine-tuning cosmological parameters (the "graceful exit" problem) [^ ■ In Brans-Dicke theory, the simplest of scalar- 
tensor theories, gravitation is described by a metric and a scalar field <f> coupled to both matter and space time 
geometry which obeys a wave equation with a source term determined by the matter distribution. Gravitational 
collapse in this theory has been studied in some of the recent literature. In this paper we study the formation of 
naked singularity as the end state of gravitational collapse of a matter fluid obeying the barotropic equation of state, 
p = wp in the context of Brans-Dicke (BD) theory. Being motivated by this concept, we investigate the violation of 
cosmic censorship conjecture in the collapse procedure of such a matter fluid. 



II. BASIC EQUATIONS 

In the context of Brans-Dicke theory with a self interacting potential and a matter field, the action is given by 

s ^hJ d * x ^ - ^v Q $v Q $ - y($)) + s m , (i) 

where the constant u> is the Brans-Dicke parameter and <I> is the Brans-Dicke scalar field, which is related to the 
gravitational "constant" by G — ^. Extremizing the action yields the field equations(we set k = 8ttG = c = 1 in the 
rest of this paper) 

Gfiv = T^ v , (2) 

where the effective stress-energy tensor is 

T^ = Ut^+1%), (3) 
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with 



and 



Tt = |(v M *v„$ - ig^V^V^) + (V M V,$ - 9llv u<s>) - \ gia ,V(*), (4) 



r™ = diag {pm,Pm,Pm,Pm) , (5) 

being the stress-energy tensors of scalar field and matter fluid, respectively. Here we look at the BD scalar field as 
matter fields originating from geometry. Variation of the action with respect to <E> gives 

™=^ + ^f*^-2V(*)V (6) 
2w + 3 2uj + 3 V d$ 7 

where T m stands for the trace of T™ and the subscript "m" refers to the matter fields (fields other that <£>). In 
general relativity the interior solution of a collapsing star is given by the Tolman-Bondi solution. Here we present 
a transparent consideration of a homogenous collapsing star with BD scalar field $ = $(t) and effective potential 
V = V($>). With this consideration, the Tolman-Bondi solution converts to a Friedmann-Robertson-Walker(FRW) 
metric. The interior metric for marginally bound (k = 0) case is given by 

ds 2 = -dr 2 + a 2 ( T )(dr 2 + rW), (7) 

where r is the proper time of free falling observer whose geodesic trajectories arc distinguished by the comoving radial 
coordinate r and dil 2 being the standard line element on the unit two sphere. Since the presence of matter acting as 
a "seed" field and origin of spherical symmetry, prompts the collapse of BD scalar field, we have considered perfect 
fluid models with barotropic equation of state being given by 

p m = wp m . (8) 

Using the continuity equation for the matter and Eq. ((5]) one gets the following relations between p m , p m and the 
scale factor as follows 

p m =P om a- 3(1+w ^, p m = wp orn a-^ 1+w \ (9) 

where p Qm = p m (a = 1), is the initial value of the energy density of matter on the collapsing shell. One then has the 
following equations for the effective stress-energy tensor 

p,.,=^-»>4('*+'™>=5('»+fl-4*+^)- <"» 

and 

Pc-n = T r ^ = T° r = T* (-> = i (p # + Pm ) = 1 L + \ % + * + 2^4 - (11) 

with all other off-diagonal terms being zero and the radial and tangential profiles of pressure are equal due to the 
homogeneity and isotropy. The interior solution of Einstein's equation for the line element (|7j! takes the form 

M' M 



R 2 = £ (13) 

The quantity M. arises as a free function from the integration of Einstein's equation and can be interpreted physically 
as the total mass of the collapsing cloud within a coordinate radius r with M > 0, and R(t, r) = ro(r) being the area 
radius for the shell labeled by the comoving coordinate r. From Eq. (|11[) we can solve for the mass function as 

M = ^(p,+ Pm ). (14) 
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Using Eqs. ([12]) and (fT"3j) we arrive at a relation between d and the effective energy density as follows 

a 2 

a 2 = +Pm), (15) 

since we are concerned with a continual collapse, the time variation of the scale factor should be negative (d < 0). 
This implies that the area radius of the shell for constant value of r decreases monotonically. For physical reasons, 
it is assumed that the energy density is non-negative every where. The singularity arising from continual collapse is 
given by a = 0, in other words when the scale factor and physical area radius of all the collapsing shells vanish, then 
the collapsing cloud has reached a singularity. A point at which the energy density blows up, the Kretschmann scalar 
K, = R ahcd Rabcd diverges and the normal differentiability and manifold structures break down. 



III. THE SOLUTION 



We would like to construct and investigate a class of collapse solutions for Brans-Dicke scalar filed with non-zero 
potential considering the matter fields, where the trapping of light is avoided till the singularity formation, thereby 
allowing the singularity to be visible to outside observers. In order to reach this purpose we consider a class of collapse 
models where near the singularity, the divergence of energy density of BD scalar field is given by the following ansatz 



P.i 



(16) 



where n is a positive constant and the scale factor, a(r), goes to zero in the limit of approach to the singularity which 
causes p^ to diverge. Using the above equation and Eq. (|15p. one may easily obtain the following relation for a as 



1 

6$ 



(2 - n)a 1-n - (1 + 3w)p gm a- (2+3w) 



6$2 



,-(1+3™) 



(17) 



where <& ja = <&/d. Now by substituting for p {a{{) and p (cff) into Eqs. (pTJj) and (fTTj) together with the use of Eqs. (|13l) . 
(IT5|) and (fTTj) we arrive at the following differential equation as 



$2 



n a 
+ 3^T 



a x - n + o a ^ (2+3w) 



5 + 3w 



3$ 2 



where we have used 



+ p a-( 1+3w *> 



6 



= 0, 



$2 

,a 

$3 



2w + 1 
6 



a 2 - ,l + p a-( 1+3w) 

i rim 



(18) 



$ = a$ a + d 2 $ aa . 



(19) 



In the following we solve Eq. (|18p for barotropic equation of state p m — wp m , where w = [0, — |, — |, ^] correspond 
to dust, cosmic string, domain wall, and radiation, respectively. Taking the following ansatz 

$(r)=a a (r), (20) 

where a satisfies the following equation 

-a 2 (l + p om )(3 + 2uj) + a(n + 4 + p am (7 + 3w)) + 2n = 0, (21) 
and by setting n = 3(1 + w), we derive an expression for the BD scalar field as a function of scale factor. 



A. Dust (w = 0) 

For such a case of pressure-less matter the parameter a will take the following values as 

12 

7(l+P om ) + V / (l+' 9 o m )( 121 + 49 '' n l + 4 8^) 

a={ . (22) 



-7(l+P oTO ) + V(l+P om )(121+49p om +48o;) 

For p m > and — § < u < oo the first value for a is negative and the second one is positive. Thus, since the BD 
scalar field must diverge near the singularity we choose the first one. For a special case in which p m = 1 and w = — 1 
(string effective action), a = —0.405. 
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B. Cosmic Strings (it) = — |) 

Cosmic strings are the consequence of 1-dimensional (spatially) topological defects in various fields. These topolog- 
ical defects are related to solitonic solutions of the classical equations for the scalar (and gauge) fields which for the 
case of a complex scalar field, cosmic strings can be formed|3(|. Such a kind of matter can be regarded as a perfect 
fluid obeying the barotropic equation of state, p m = — l/3p m . We are interested here to find the behavior of BD scalar 
field as a function of scale factor in gravitational collapse of this kind of fluid. The corresponding values of a for this 
case are 



{3(l+P om ) + v / (l+P om )(21+9p om +8o;) 
(23) 
-3(l+P om ) + VU+P,,™ )(21+9p Qm +8u) 

It is seen that for p m > and — | < u < oo the first and second values of a are negative and positive, respectively. 
Choosing the first value and setting p Qm =1 and oj = — 1, we have a = —0.316. 



C. Domain Walls (w 



Domain walls are Generally topological solitons which occur whenever a discrete symmetry is spontaneously 
broken[3lJ. These are the result of 2-dimensional topological defects in different scalar or gauge fields. As for 
the case of cosmic strings, domain walls can be regarded as a perfect fluid obeying the barotropic equation of state, 
p m = — 2/3p m . We here again try to find the behavior of BD scalar field as a function of scale factor in gravitational 
collapse of such kind of matter fluid. In order to reach this purpose we consider the following values obtained for a as 



-HT-+P om ) + y/(T-+P om )(^+25p ()m +Wu) 

(24) 

' 5(l+P nm ) + V( 1 +P om )( 4 9+ 2 5P om +16") 

As it is seen the second one is always negative for p Qm > and — | < u> < oo causing the BD scalar field to diverge 
in the vicinity of singularity. For uj = —1 and p Qm = 2 one has a = —0.129. 

D. Radiation (w = |) 

Finally for this case of matter, we find the following values for a as 

4 



(25) 
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-^+P om ) + y / ^+P om )(5+2p om +2uj) 

Again the first value is always negative for p m > and — | < ui < oo. Setting p m = 2 and oj = — 1 we have 
a = -0.32. 



IV. TIME BEHAVIOR OF THE SCALE FACTOR 



One would like to study the time-dependence behavior of the scale factor during the collapse procedure, considering 
matter. If at time r = r* (or equivalently for some a = a*) the energy density of the BD scalar field starts growing 
as a - ™, then by integrating Eq. (|15[) in the vicinity of the singularity with respect to time one gets the time behavior 
of the scale factor as 

a(r)= /a**(°+3U+™)) _ I sj^f^ (a + 3(1 + w)){r-rM \ (26) 



and the corresponding singular epoch as 



= 2x/3 

s y /T+^( a + 3(1 + w)Y 1 ' 

where the time t s corresponds to a vanishing scale factor. Thus the collapse reaches the singularity in a finite proper 
time. This result for the scale factor completes the interior solution within the collapsing cloud, providing us with 
the required construction. 

V. CONDITIONS ON RADIAL NULL GEODESIC EXPANSION 

Consider a congruence of outgoing radial null geodesies having the tangent vector (£ r , £ r , 0, 0), where £ T = 
dt/dX and £ r = dr/dX and A is an afhne parameter along the geodesies. In terms of these two vector fields the 
geodesic equation can be written as 



d? _ 2a cter 



dX 

and 



d^_ 
~dX 



eC, (28) 
a 



ad(C) ■ (29) 



The geodesic expansion parameter is given by 

<-> v ; r dtf • r;,r. (30) 



which gives 



dr dr 
In order to compute the sum 



we proceed by noting that 



and similarly, 



dr dr 



dC _ di T dr dC dr 
dX dr dX dr dX' 



dC _ dC dr d£ r dr 
dX dr dX dr dX 



(32) 



(33) 



(34) 



Dividing the first of these two relations by dr/dX and the second by dr/dX, and after adding the resulted equations 
one gets 
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where we have used Eqs. (|28|) and (J29J) , and the fact that for outgoing radial null geodesies the relation between £' 
and £ r is given by 

F = £ = " M ' (36 » 

Substituting Eq. ([55]) into Eq. (|3"TT) and after a simple calculation we arrive at the desired expression for 0: 



e = ;U-\/^ )■ ( 37 ) 

In order to determine the visibility, or otherwise, of the singularity, one needs to analyze the behavior of non-spacelike 
curves in the vicinity of the singularity and the causal structure of the trapped surfaces. These surfaces are closed 
oricntable smooth two-dimensional space-like surfaces such that both families of ingoing and outgoing null geodesies 
orthogonal to them necessarily converge (32[. The singularity will be called naked if there exists a family of future 
directed non-spacelike geodesies, reaching faraway observers in space-time and terminating at the singularity in the 
past. The existence of such curves implies that either photons or time-like particles can be emitted from singularity. 
So if the null geodesies terminate at the singularity in the past with a definite tangent, then at the singularity we 
have > 0. If such family of curves do not exist and the event horizon forms earlier than the singularity covering it, 
a blackholc is formed. The boundary of the trapped surface region in the space-time is called apparent horizon where 
in spherically symmetric space-time is given by 

//"'//,//,, = 0. (38) 

Therefore at the boundary of the trapped surface the vector is null. Using Eq. (JT]), the above equation can be 
written as 

- R 2 + a- 2 R' 2 = 0, (39) 

which leads to M. = R. Here use has been made of Eq. (|13[) . The space-time region where the mass function A4 
satisfies A4 < R is not trapped, while M. > R describes a trapped region 0, H^, HH . 

Let us now study the relation that O bears with the formation or otherwise of a naked singularity in spherical 
collapse. Calculating A4/R in the general case which is considered for energy density and by using Eq. (|13p . we 
have 

r a 2 - n + p nm a-^ +3w A. (40) 



R 3$ 

We shall employ the above equation and Eq. (|37p to examine the nakedness of the singularity as the collapse end state 
for the four cases of matter field considered in section III. We show that physically, the formation of a black hole or 
a naked singularity as the final state for the dynamical evolution is governed by the rate of collapse scenario and the 
presence of BD scalar field. It is seen that the cosmic censorship conjecture is violated for w = — | and w = — |. The 
weak energy condition which states that the energy density as measured by any local observer must be non-negative 
can be written for any timclikc vector as follows 

T^V»V U > 0, (41) 

whereby one gets the following conditions for the effective energy density (p (off) > 0) 

(l + p om >- (Q+3(1+,u)) >0, (42) 
and the sum of effective energy density and pressure (p (off) + p (oft) > 0) as 

w = 0, a < -> [(1 + p Qm )(l - f )] a~( Q + 3 ) > 0, 

w = -i, a < -> [(1 + p om )(^)] a~ (Q+2) > 0, 

(43) 

w = -§, a < [(1 + p om )(^)] a-(« +1 ) > 0, 
w = i, a < -> [(1 + p om )(^)] a- (Q+4) > 0. 
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It is seen that these two conditions are satisfied by all cases of w and p g , n > and — | < oj < oo considered above. 
At the initial epoch (a = 1) there should not be any trapping of light. Assuming r = is the boundary of the 
collapsing ball, then at the initial epoch the ratio Ai/R is less than unity for a suitable initial value of p gTa standing 
for all values of w. This fact is in accordance with the regularity condition stating that the gravitational collapse 
must initiate from regular and physically reasonable initial conditions. The time at which the physical area radius of 
the collapsing cloud becomes zero, denotes a shell- focusing singularity which lies on the curve R(r s ,r) = where t s 
being the singular epoch given by Eq. (|27[) . For the case of homogeneous-density collapse the resulting singularity 
may lay on the curves i?(r s ,0) = or R(r s ,r =/= 0) = 0, which corresponds to a central or non-central singularity, 
respectively We consider first the simpler case of non-central singularity and investigate the failure of formation of 
apparent horizon in collapse scenario for different values of w. 



A. Dust (w = 0) 

We are now in a position to study the effect of BD scalar field on the formation or otherwise of the apparent horizon 
as the dynamical procedure of collapse scenario evolves ( we set u> = —1 in the rest of this paper). We begin by Eq. 
(fT4)) which for w = can be written as 

^ = y(l + /V>-< 1+ «>. (44) 

The initial energy density of matter must be positive due to the regularity conditions, then for p gTn > and oj = — 1, 
the first value of a in Eq. (f22|) implies that |a| < 1. From Eq. (|44|) it is seen that the ratio M./R grows and the 
expansion parameter, Eq. (|37|) . tends to negative infinity. Thus there exist no radial null geodesies emerging from 
the singularity. Strictly speaking the singularity occurred here is necessarily covered and a black hole is formed as 
the collapse end state. 



B. Cosmic Strings (w — — |) 

For this case Eq. (fbf)) and the time variation of the mass function take the following form as 

(45) 

M = -(l + p om )^{a-l)aa- a . 

As it is seen from Eq. ([23]) . since the first value of a is always less than zero the ratio M./R stays finite till the 
singular epoch and causes the expansion parameter to be positive up to the singularity, and if no trapped surfaces 
exist initially then no ones would form until the epoch a(r) = which is consistent with the fact that there exist 
families of outgoing radial null geodesies emerging from the singularity. One can take the positive value of a in Eq. 
(j23|) . but for this case the BD scalar field get vanished as the scale factor tends to zero. Also the weak energy condition 
may be violated. In addition to, for such a case the ratio M./ R grows at a vanishing scale factor causing the expansion 
parameter tends to negative infinity which means that the singularity is covered and no radial geodesies can emerge 
from it. From the second equation in Eq. (|45j) it can be seen that the time derivative of the mass function for a < 0, 
is negative (note that a < 0) which means that the mass function contained in the collapsing shell with that radius 
keeps decreasing. In other words there exists an outward energy flux during the collapse scenario. Since no trapped 
surfaces form up to the singularity, the outward energy flux would be observable. 



C. Domain Walls (w = -§) 

For this case one may rewrite Eq. (1141) and time derivative of the mass function as 

R 3^ i + P m > ' 

(46) 

M = -(l + p Qm )^(a-2)aa 1 - a . 

from the first equation one may easily see that at initial epoch (a = 1), the regularity condition is satisfied. Since the 
second value of a in Eq. (j2~4")l is always negative, the ratio of mass function to area radius of the collapsing shell is less 
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than unity during the collapse procedure denoting that the expansion parameter being positive up to the singularity. 
In this case the collapse evolution to a naked singularity takes place, where the trapped surfaces do not form early 
enough or are avoided in the spacetimc. For first value of Eq. (f24|) . a > 5 which causes the ratio Ai/R tends to 
infinity as the scale factor vanishes and goes to negative infinity, thus trapped surfaces do form in the spacetime 
which prevent the null geodesies to emerge from the singularity. Such a situation ends in a black hole as the final 
fate of the collapse scenario. But such a value of a is not allowed since it violates the weak energy condition. From 
the second equation in Eq. (|46l) . it is obvious that for negative value of a and d < 0, the time derivative of the mass 
function is negative stating that the mass contained in collapsing ball reduces as the time advances. 



D. Radiation (w — 



In this case Eq. ([14"]) can be written as 



M 
R 



(47) 



From Eq. (|25[) one can easily see that the first value for p o , n and ui — — 1 is always negative and \a\ < 4. Thus, 
in such a situation the ratio M. j R tends to infinity as the singularity is approached. Thus the expansion parameter 
behaves just as the dust case, and the final singularity is necessarily covered within an event horizon of gravity. 

The central singularity occurring at R = 0,r = is naked if there exist outgoing non-spacelike geodesies reaching 
faraway observers and terminating in the past at the singularity. In order to investigate the nakedness of this kind 
of singularity we proceed by introducing a new variable x = r s 1 and 8 > 1 is defined such that iJ'/r" 5-1 is a unique 
finite quantity in the limit r — > 0. Then we have the following equation 



dR 
dx 



1 



a-l \ fa 



R' 



By virtue of Eqs. (fT3"|) and (|3"o]l the above equation leads to 



dR 

dx 



Rf 



ax 



1- 



(48) 



(49) 



It is clear that i? = 0,a; = 0isa singular point of Eq. (|4"9"|) . If there are outgoing radial null geodesies terminating in 
the past at the singularity with a definite tangent, then at the singularity we have 4^ > 0. For w = — h, and w = — | 
with a being negative the quantity A4/R < 1 throughout the collapse procedure, so the term being in the second 
bracket is positive and 4^ > as the singularity is approached indicating that the singularity is visible to outside 
observers and the inverse result holds for w — and w = g. 



VI. NAKEDNESS OF THE SINGULARITY 



The continual gravitational collapse of a matter cloud culminates in either a black hole or a naked singularity 
where in the former an event horizon develops earlier than the formation of the singularity. Thus the regions of 
extreme physical conditions such as densities and curvatures are hidden from the outside observers. If such horizons 
are delayed or failed to develop during the collapse procedure, as governed by the internal dynamics of the collapsing 
cloud, then the scenario where the ultra-strong gravity regions become visible to external observers occurs and a 
visible naked singularity forms. In such a case where no black hole forms, the field collapses for a while and then 
disperses. Therefore as viewed by a central observer, the scalar invariants namely the Kretschmann scalar should 
grow near the singularity, gain some maximum value and then approach to zero at late times [35J. Since the absence 
of an apparent horizon does not necessarily implies the absence of an event horizon, we examine the nakedness of the 
singularity in spherically symmetric collapse of a fluid by considering the behavior of the Kretschmann invariant with 
respect to time. For the line element (7) this quantity is given by 

JC = R abcd R abcd = H [ a V + d 4] _ (50) 
a 4 

By the virtue of Eqs. (|26p and (|2T1) for w = — | the above quantity can be written as 

3.322 

^ cs (l-0.6r)4> 
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and for w = — | as 

Kdw = (l-O.tr)*- (52) 

Note that both these two cases satisfy the condition on expansion parameter stating that this quantity must be 
positive up to the singularity. Fig. 1 and Fig. 2 show the behavior of the Krctschmann scalar as a function of proper 
time, t. It is seen that both K, cs and K-dw diverge at r = 5/3 and r = 5/2, respectively. They then converge to zero 
at late times signaling the failure of formation of the event horizon. 

Let us now consider the geometry of the exterior spacetime. In order to fully complete the spacetime model, one 
needs to match the interior spacetime of the dynamical collapse to a suitable exterior geometry. The Schwarzschild 
solution is a useful model describing the spacetime outside the sun and stars. However this model may no longer 
be suitable to describe the exterior geometry of any realistic star, because the spacetime outside such a star may be 
filled with radiated energy from the star in the form of electromagnetic radiation. The Schwarzschild model does 
not describe this as it corresponds to an empty spacetime given by T a b — 0. The spacetime outside a spherically 
symmetric star being surrounded by a radiation emitted from the star is described by the Vaidya metric [36[ which 
can be given in the form 



ds'i ut = - I 1 - 2M ^ r »' u ) ^ du 2 _ 2dudr u + rl<m 2 , (53) 

where u, being the retarded null coordinate, r u and M(r u ,u) are the Vaidya radius and Vaidya mass, respectively. 
Following the work of [33[ we use the Isreal-Darmois junction conditions to match the interior spacetime described by 
Eq. ([7]) to a Vaidya exterior geometry at the boundary hypersurface E given by r = rj,. The spacetime metric just 
inside E is given by 

ds 2 n = -dr 2 + a 2 (r) [dr 2 + r b W] (54) 
Matching the area radius of the collapsing shell at the boundary, one gets the following equation 

r u (u) = r b a(r), (55) 
whereby on the hypersurface E, the interior and exterior metrics can be written as 

dsl in = -dr 2 +a 2 (T)r 2 dn 2 , (56) 

and 

„2 _ (, 2M(r„,«) dr u \ , , 2 



dsl out = ( 1 - ; + 2^ ) du 2 + r 2 u dQ 2 . (57) 



Matching the induced metric on £ one gets 

fdu\ 
[d^J 



( du\ 1 , . , , , „ 

hr = r. (Tuh = na(r). (58) 



S {-y 2M(r u ,u) , r,dr l 



In order to match the extrinsic curvature for interior and exterior spacetimes, one has to fine the unit normal vector 
field to the hypersurface E. In this step, we proceed by noting that any spacetime metric can be written locally in 
the form 

ds 2 = - (N 2 - NiN 1 ) dr 2 - 2N l dx l dr + h^dx'dx 3 , (59) 

where N, N l , and hij are the lapse function, shift vector, and induced metric, respectively and i, j are three-dimensional 
indices run in {1,2,3}. The contravariant and covariant components of the unit normal vector field are given by 

n a = ± (S a - N a ) , n a = -NS° a . (60) 

Comparing Eqs. (|59p and ([7)1. one finds the contravariant components of the normal to the hypersurface E for the 
interior metric as 

<„ = [0,a(r)- 1 ,0,0]. (61) 
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Upon using a similar approach, one finds the first non-vanishing contravariant component of the normal to £ for the 
exterior metric as 

n u = r . (62) 

^ _ 2M(ry, M ) 9 dr, ' 
r u du 

In order to compute the second non-vanishing contravariant component of the normal vector field we proceed by 
having recourse the normalization relation holding for n a as 

n u n u + n Tu n ru = 1. (63) 

Benefiting from the property of the metric tensor in raising and lowering indices, we have the following relations 

_ 2M{r u ,u) 

-n r "; n ru =-n u . (64) 



_ 2M(r u ,u) n dr u 
r u ^ du 



Substituting the above equations and Eq. ([6"2")l into Eq. (|6"3")l and after a simple calculation we arrive at the desired 
expression for n r ™ as 

_ 2M(r u ,u) dr u 

n " = — . (65) 

j _ 2M(r„,n) r,dr u ' 

The extrinsic curvature of the hypcrsurface £ is defined as the Lie derivative of the metric tensor with respect to the 
normal vector field being given by the following relation 

K ab = ^ [gab, c n c + 9cbn c a + ga C n c b } ■ (66) 

Since the matching is for the second fundamental form, K a t, there exists no surface stress energy or surface tension 
at the boundary |37|. The nonzero 9 components of the extrinsic curvature are 

^ _ 2M(r u ,u) , dr^ 

K& = r b a(r) , K° 9 f = r u ^ . (67) 

^ _ 2M(r„u) , 2 ^r„ 
r u du 

Setting [Kgg — Kgg 1 ] = on the hypersurface S, and by using Eqs. (fl"3"|) and (|55)l. one gets the following relation 
between mass function and Vaidya mass on the boundary as 

M(T,r b ) = 2M(r u ,u). (68) 

From the above equation and Eq. ([T4"]) one can see that the BD scalar field affects on the behavior of the Vaidya mass 
in the collapse scenario. In order to find another relation expressing the rate of change of Vaidya mass with respect to 
r u one has to match the r component of the extrinsic curvature on the hypersurface E. Having set [K™ — K°™ 1 ] = 0, 
one gets 

M {r u , u) ru = — + r 2 b aa. (69) 

The occurrence of a naked singularity as the final outcome of a collapse procedure, depends on the non-existence of 
trapped surfaces till the formation of the singularity, which corresponds to the existence of families of non-spacelike 
trajectories reaching faraway observers and terminating in the past at the singularity. In order to show this, we begin 
by Eq. (|58p and after using Eqs. ([67)) . and (f68|) we arrive at the following relation 

du\ l-r b a 

dr)^ 1 M(r,r b ) ^ 

It can be easily checked that if one imposes the null condition on the Vaidya metric, the result is the same as Eq. [70] 
What is meant by this, is that null geodesies can come out from the singularity and reach to faraway observers before 
it evaporates into the free space. In other words the formation of trapped surfaces in spacetimc is avoided and a 
naked singularity can be produced. 
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VII. BEHAVIOR OF THE BD SCALAR FIELD AND THE EFFECTIVE POTENTIAL 



In the following section we wish to study how the effective potential behaves as the scalar field varies. For this 
purpose we start by Eq. (jTU)) and consider the four cases of matter field discussed in section III. Together with the 
use of Eqs. (fT5|) and ([2TH) . one may easily find the following expression for the potential as 



7(0)=^"^, p= (2+^(1 + Pom )(6-c) 



(71) 



Fig. 3 shows the behavior of the BD potential with respect to 4> for different values of uu and u = — 1. 

Let us now consider the case in which the BD scalar field is a function of both r and r. Assuming that far away 
from the collapsing system, the effective energy density behaves homogeneously, we obtain a measure of the radial 
profile of the BD scalar field for each cases of w considered in Section III. We begin by Eq. © together with the use 
of Eqs. (|T5l) . (|T7|) . (|T9|) . and (|7T1) we arrive at a differential equation for $(a(r), r) as 



6 6 



where C is given by 



,-(1+3™) 



a 2 ra 2 2a; + 3 



a 2 -™+p m a- (1+3lu) 



C 



3/3(1 + 10) 



+ 2/3-p om (3 W -l). 



3^ 



(72) 



(73) 



the above equation can be more simplified and the result is as follows 



- Da 



l—a—3w 



= 0, 



(74) 



where D is a constant and is given by 



5a — a 2 aw\ , . C 

— g— --2-)(i + p,J-aJ+ 



(75) 



Here / denotes partial differentiation with respect to r. Solving the above differential equation with a suitable choice 
of value for p m one may find the solutions as functions of t and r. Figs. 4-7 show the behavior of BD scalar field in 
terms of r and r in which the constants of integration have been chosen in such a way that the scalar field blows up 
near the singularity and get vanished faraway from the collapsing system. 



VIII. CONCLUSION AND OUTLOOK 



In this work we have studied the gravitational collapse of the BD scalar field with non-zero potential in the presence 
of matter fluid. Assuming that the energy density of the BD scalar field behaves as the inverse power law of the scale 
factor near the singularity, we presented a class of solutions in Brans-Dicke theory in which the naked singularity 
can be created being accompanied by the violation of the cosmic censorship conjecture. In section III, we found 
the behavior of BD scalar field as a function of scale factor for the four cases of matter fluid. In section IV, the 
general expressions of time behavior of the scale factor and singular epoch has been achieved. Having examined the 
ansatz taken for divergence of the energy density of BD scalar field, Eq. (|16p . together with the use of Eq. ((9]) for 
energy density of matter fluid and by using the concept of expansion parameter, we have shown in section V that the 
presence of BD scalar field due to the Eq. (fT^|) can affect the formation or otherwise of the trapped surfaces and only 
for the two cases, w = — 4 and w = — |, formation of the apparent horizon can be failed and a naked singularity may 
be generated as the final fate of collapse procedure. But since the absence of an apparent horizon does not necessarily 
implies the absence of an event horizon, we have computed the Kretschmann scalar in section VI and the result has 
been plotted in Fig. 1 and Fig. 2 for w = — -| and w = — |, respectively. It is seen that this quantity diverges 
at singular epoch, and then vanishes at late times, a behavior which can be interpreted as the absence of an event 
horizon and formation of a naked singularity. Also following the work of [33| we have shown at the end of this section 
that the Vaidya geodesic emerging from the singularity before it evaporates into free space is null. 
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Beside our work which have only treated exact solutions to gravitational collapse in Brans-Dicke theory, one may 
find numerical solutions to such a collapse scenario in the literature In HEHII, tne aut hor has developed a new nu- 
merical code that solves the gravitational field equations coupled to the matter for evolution of a spherically symmetric 
configuration of noninteracting particles in Brans-Dicke theory. Using this code, he has shown that Oppenheimer- 
Snyder collapse in this theory results in black holes rather than naked singularities, at least for |3 + 2lj\ > 3, which 
are identical to those of general relativity in final equilibrium, but are quite different from those of general relativity 
during dynamical evolution in which they radiate mass. The reason for this behavior is due to the violation of the null 
energy condition even in vacuum spacetimes with positive values of lo, passing of the apparent horizon of a black hole 
outside the event horizon, and decreasing the surface area of the event horizon over time. This numerical code enables 
one to decide on a number of long-standing theoretical questions about collapse in Brans-Dicke theory of gravitation. 
Also the gravitational collapse of a scalar field with other characteristics and couplings has been discussed in some 
literature. In [3i| , the collapse of a self-similar scalar field has been studied, and it has shown that there exists two 
classes of solutions which one of them consists of a nonsingular origin in which the scalar field collapses and disperses 
again. There is a singularity at one point of these solutions which is not observable at a finite radius. The second class 
of solutions contains both black holes and naked singularities with a critical behavior interpolating between these two 
extremes. Numerical study of spherically symmetric collapse of a masslcss scalar field has presented in (ioj . where it 
is shown that the masses of black holes which form satisfy a power law Mbh oc \p — p*\ J ■ Where p is a parameter 
which characterizes the strength of initial condition, p* is the threshold value and 7 ~ 0.37 is a universal exponent. 
Also the collapse of a massless scalar field in Brans-Dicke theory has studied both analytically and numerically in 
fill ] and it is shown that for oj > — |, a continuous self-similarity continues and that the critical exponent 7 depends 



on u>. In [42(, gravitational collapse of a self-interacting (massive) scalar field has been studied both analytically and 



numerically on a Reissner-Nordsrom background and finally in [43j , one may find some examples of naked singularity 
formation in gravitational collapse of a scalar field. 
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FIG. 2: The behavior of Kretschmann scalar (in units of s 4 ) as a function of proper time for u) = —1 and w 
density, scale factor, and proper time we have adopted the values p„ = 2, a* = 1, and r* = 0, respectively. 
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FIG. 3: Behavior of potential as a function of BD scalar field for lj = —1 and different values of p and 
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FIG. 4: Behavior of the BD scalar field with respect to r and r for u) = —1, w 
proper time we have adopted the values p m — 1, a* = 1 and t* = 0, respectively. 
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FIG. 5: Behavior of the BD scalar field with respect to r and r for u = —1, w = 0. For the initial energy density, scale factor and proper 
time we have adopted the values p m = 1, a* = 1 and t* = 0, respectively. 




FIG. 6: Behavior of the BD scalar field with respect to r and r for uj = —1, w = — j|. For the initial energy density, scale factor and 
proper time we have adopted the values p m = 2, a* = 1 and t* = 0, respectively. 
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FIG. 7: Behavior of the BD scalar field with respect to r and r for ui = —1, w = i. For the initial energy density, scale factor and proper 
time we have adopted the values p = 2, a* = 1 and t* = 0, respectively. 



